The Scaling of the Anomalous Hall Effect in the Insulating Regime 



Xiong- Jun Liu, Xin Liu^ and Jairo Sinova^^ ^ 

^Department of Physics, Texas A&M University, College Station, Texas 77843-4242, USA 
^Institute of Physics ASCR, Cukrovarnick 10, 162 53 Praha 6, Czech Republic 

(Dated: January 19, 2013) 

We develop a theoretical approach to study the scaling ol anomalous Hall effect (AHE) in the 
insulating regime, which is observed to be oc cr^jf" ^^ '^^ in experiments over a large range of 

materials. This scaling is qualitatively different from the ones observed in metals. Basing our theory 
on the phonon-assisted hopping mechanism and percolation theory, we derive a general formula 
for the anomalous Hall conductivity, and show that it scales with the longitudinal conductivity 
as a'^^ ~ a],x with 7 predicted to be 1.38 < 7 < 1.76, quantitatively in agreement with the 
experimental observations. Our result provides a clearer understanding of the AHE in the insulating 
regime and completes the scaling phase diagram of the AHE. 
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The anomalous Hall effect (AHE) is a central topic in 
the study of Ferromagnetic materials [T] . It exhibits the 
empirical relation p,j.y = RqBz + RsM^ between the total 
Hall resistivity and the magnetization and external 
magnetic field B^. Here Rq and Rs are respectively the 
ordinary and anomalous Hall coefficients. When trans- 
formed to an anomalous Hall conductivity (AHC), , 
three regimes are observed with respect to its dependence 
on the diagonal conductivity, axx- In the metallic regime 
the AHE is observed to be linearly proportional 

to axx for the highest metallic systems {axx > lO^fJ^^ 
cm~^) and roughly constant for the rest of the metallic 
regime. This dependence indicates the different domi- 
nant mechanisms in ferromagnetic metals. These are un- 
derstood to be the skew scattering, side jump scattering, 
and intrinsic deflection mechanisms. The intrinsic con- 
tribution is induced by a momentum-space Berry phase 
linked to the electronic structure of the multi-band SO 
coupled system [Tl[2]- The side jump scattering mecha- 
nism gives the same scaling relation as the intrinsic con- 
tribution, i.e. ax^~^^ oc cr°^, and the skew scattering is 
linear in the longitudinal conductivity, a^^^'^^ cx axx- 
While these mechanisms are now better understood, the 
maximum scaling exponent of the AHC cannot exceed 
unity in the metallic regime [1]. 

On the other hand, experiments in the insulating 
regime exhibit an unexpected scaling relation of the AHC 
ax^ (X a^.'x'^""^''^^ , which remains unexplained and a ma- 
jor challenge in understanding fully the phase diagram 
of the AHE [IHSj . The available microscopic theories of 
metals fail in this regime since the condition kpl ^ 1 is 
no longer satisfied for disordered insulators [TJ [T3] . The 
few previous studies of the AHE in the insulating regime 
focused on manganites and Gai_a;Mn2;As; while the man- 
ganites do not exhibit this scaling, the studies on insu- 
lating Gai_2,Mn2,As did not show this scaling [TMT5] . 

In this Letter we study the scaling of the AHE in the in- 
sulating strongly disordered amorphous regime, where at 
low temperatures charge transport results from phonon- 



assisted hopping between impurity localized states [El 
f20| . We calculate the upper and lower limits of the 
AHC, and show it scales with axx as a^^ ~ a^^ with 
7 predicted to be 1.38 < 7 < 1.76, in agreement with 
the experimental observations. This scaling remains the 
same regardless of whether the hopping process is Mott- 
variable-range-hopping or influenced by interactions, i.e. 
Efros-Shkolvskii (ES) regime. 




FIG. 1: AHE in the insulating regime. In this regime charge 
transport occurs via hopping between impurity sites. 

To capture the Hall effect one requires the hopping 
process between impurity sites (Fig. [I]) to break the time- 
reversal (TR) symmetry. The two-site direct hopping 
preserves TR symmetry, and therefore more than two 
sites must be considered. The hopping through three 
sites, as depicted in Fig. [2] is the minimum requirement 
to model theoretically the ordinary Hall effect (OHE) 
|21) . The total hopping amplitude is obtained by adding 
the direct and indirect (through the intermediate fc-site) 
hopping terms from i to j sites. The two hopping paths 
give rise to an interference term for the transition rate 
which breaks TR symmetry and is responsible for the 
Hall current in the hopping regime. For the OHE, the 
interference is a reflection of the Aharonov-Bohm phase, 
and for the AHE it reflects the Berry phase due to SO 
coupling. Furthermore, the dominant contribution to the 
Hall transport will be given by the one- and two-real- 



2 



phonon processes through triads (Fig. [2| pT] . 

Our theory is based on a minimal tight-binding Hamil- 
tonian. With the particle-phonon couphng considered, 
the total Hamiltonian H = Hp + He + Hph , with 

iaX 

Hph = y^^u!\b\bx. 
A 

Here Hp describes localized states, He gives the particle- 
phonon coupling with 77 the coupling constant, Hp^ is the 
phonon Hamiltonian, a is the local on-site total angular 
momentum index, and is the energy measured from 
the fermi level. Here we consider that the magnetization 
is saturated and thus assume M = Mcz- The hopping 
matrix tij is generally ofF-diagonal due to SO coupling 
(see Supplementary Information (SI)). The localization 
regime has the condition l^iaj^l <C |ei — ej| in average. 
The specific form of the relevant parameters {tij, M, spin 
operator Tap) are material dependent and do not affect 
the scaling relation between and a^x- 
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FIG. 2: (Color online) The hopping processes through triads 
with up to two real phonons absorbed or emitted. (Top) Typ- 
ical diagrams of the two-phonon direct and indirect hopping 
processes. (Bottom) One-phonon direct process and typical 
three-phonon (one real phonon) indirect hopping processes. 

Considering the dominant contributions to the longi- 
tudinal and Hall transports, we obtain the charge current 
between i and j sites in a single triad with applied volt- 
ages [17]: lij = GijVij+Q^jk{V^k + Vjk), with the di- 
rect conductance and Qijk responsible for Hall transport. 
The formula of Uj gives the microscopic conductances 
in any single triad (see SI). To evaluate the macroscopic 
AHC, we need to properly average it over all triads in 
the random system. This is achieved with the aid of per- 
colation theory, a fundamental tool to understand the 
hopping transport. 

We first map the random impurity system to a ran- 
dom resistor network by introducing the connectivity be- 
tween impurity sites with the help of a cut-off Gc(T). 



When the conductance between two impurity sites sat- 
isfies Gij > Gc, we consider the i,j sites are connected 
with a finite resistor — 1/Gy . Otherwise, they are 
treated as disconnected, i.e. Gij — >■ 0. The Hall effect 
will be treated as a perturbation to the obtained resis- 
tor network. The cut-off Gc should be properly chosen 
so that the long-range critical percolation paths/clusters 
appear and span the whole material, and dominate the 
charge transport in the hopping regime. The macroscopic 
physical quantities will finally be obtained by averaging 
over the percolation path/cluter. 

The hopping coefficient generally has the form tia,jf3 = 

^ia jl3^~'^^^^ ^ with a^^ the localization length and Rij = 
|Ri — Rj|. The direct conductance holds the form Gij — 
Go(T)e-2'^^'^-5/3(l^-l+l^./3l+l<^--^j/3l), and then the cut- 
off can be introduced by Gc ~ Gqc^^^'''^'^^ Here /3fc 
is a decreasing function of T, indicating the material in 
the insulating regime. The number of impurity sites con- 
nected to a specific site i with energy can be calcu- 
lated by n{€i,^c) = J dej J d^R^jpicj, Ri)Q{Gij ~ Gc)- 
Here 6(a;) is the step function and the DOS p{e,Ri) w 
^^,-(5(e — ei) is approximated to be spatially homo- 
geneous. The number n{ei,^c) can also be given by 
niei^^c) = Z]«-P"(<^»;Cc), with P„(ej,^c) being the proba- 
bility that the n-th smallest resistor connected to the site 
i has the resistance less than l/Gc- The function P„ reads 
Pn{f^i,^c) = („^i)! Jo"'''*'' e"^a;""Ma; 0. The percolation 
path/cluster appears when the average connections per 
impurity site n ~ {n{ei))c reaches the critical value fie, 
where the definition of (...)c is given in Eq. ([6|. Suppose 
a physical quantity F{ei, ...,6™; ri, r„i) being a m-site 
function, requiring the i-th site to have at least rji sites 
connected to it. The averaging of F{e; f) reads 

{F{e;r))c J dei... J de,n J d^ri2... J d^r,n_i^jn 

m 

X n ^').(<^»)-^(^l' ^1' ^m), (1) 

1=1 

where Afp is a normalization factor and the probabil- 
ity function Vmi^i) = P[^i)Y.k>r,,Pk{>^i)- The term 
X]fe>i( Pki^i) entering the probability function has im- 
portant physical reason. The configuration averaging is 
not conducted over the whole impurity system, but over 
the percolation cluster which covers only portion of the 
impurity sites. Therefore the probability that an im- 
purity site belonging to the percolation cluster must be 
taken into account for probability function. Moreover, 
this probability function also distinguishes the physical 
origins of the AHC and a^x- For one has 77^ = 3, 
and for Uxx one has rji — 2. This indicates the averag- 
ing of Gxx is performed along the one dimensional (ID) 
percolation path, while for AHE which is a two dimen- 
sional (2D) effect, one shall evaluate AHC over all triads 
connected in the 2D percolation cluster. 
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Numerical solutions show the critical site connectiv- 
ity is fie = 2.6 ~ 2.7 for the appearance of a percolation 
path/cluster in three dimensional materials [MlEl]. This 
indicates the triads are sparsely distributed in the per- 
colation cluster, as shown in Fig. [3] The AHC can be 
derived by examining the transverse voltage V,^ (along 
the y-axis) induced by the applied longitudinal current 
Jq. Denote by N{x) the number of triads distributed 




FIG. 3: (Color online) Typical resistor network in the ma- 
terial. The present situation indicates V^_2 and in the 
region from x — Ax to a; -I- Aa; are zero, where no triads form. 

along the y-axis in the region around position x (Note M 
is along the z-axis, hence we assume the system in this 
direction to be uniform). The transverse voltage equals 
the summation over the voltage drops of the N{x) triads: 

Vy"ix) = 'EfJi^ y"- The average Hall vohage can 
be obtained in the limit N{x) 
we find (see SI for details) 



oo, which from Eq. Q 



AH 



2 ksT J2a0-y [i-MUa,jf3tjp,kftk-(.ia)T-^l] , 
l^i^j^k I'-u'-jfcl -^ij -^jk 



,(2) 



with L the correlation length of the network. Note the 
configuration integral given by Eq. ^ is first derived for 
the AHC in this letter. This is an essential difference from 
the former theory by Burkov et al |17j . where the con- 
figuration averaging applies to the whole system rather 
than to 2D percolation cluster. With our formalism the 
key physics that Hall currents are averaged over perco- 
lation clusters can be studied, which is a crucial step to 
understand the insulating regime of the AHE phase dia- 
gram. The above configuration integral cannot be solved 
analytically. In the following we study the upper and 
lower limits of the AHC by imposing further restrictions 
in Eq. ([2|, with which the range of the scaling relation 
between and axx can be determined. 

The lower (upper) limit of the AHC can be formulated 
by keeping only the maximum (minimum) term in the 
denominator and the minimum (maximum) term in the 



numerator. Furthermore, for simplicity we first approxi- 
mate the DOS to be constant although this approxima- 
tion is relaxed later. As a result, with further simplifica- 
tion (see SI) we find 



y j min 
max 



eH 



(0) 

max / rain 



_ I i:>max\ I ^max\ 



(3) 



— pa.(Rij+Rjh 



where (i?™r)c = e" 

g0.5/3(|ei|-H|ej|-H|£i-efc|-|e,-efc|)c 



"r)n and 



^i7fc^)c hold the same form for the calculation but the 



restrictions change to be Rij,Kjk > Kik and \ei\ > \ej\ > 



|efc|, respectively. The coefficient t 



Rik > Rik and 

(0) 

max/min 



represents the 

(0) 



maxmimum/minimum element in the matrix t]; . It is 
instructive to point out the underlying physics of the two 
limits. In the hopping regime, charge transport may pre- 
fer a short and straight path in the forward direction 
with larger resistance than a long and meandrous path 
with somewhat smaller resistance [21|TH]. This picture in- 
troduces an additional restriction complementary to the 
percolation theory for charge transport. What bonds 
in a triad play the major role for the current fiowing 
through it is determined by the optimization of the resis- 
tance magnitudes and spatial configuration of the three 
bonds. A quantitative description can be obtained by 
phcnomenologically introducing an additional probabil- 
ity factor to restrict the charge transport [H [TS] . Here 
we only need to adopt this picture to present the two 
extreme situations corresponding to {<yxy^}min/max- To 
get the upper limit we assume that for each triad of the 
percolation cluster the two bonds with smaller direct con- 
ductance dominate the charge transport, i.e. the product 
of two smallest conductances minimize the denominator, 
and take the maximum value for the numerator of Eq. 
([2]). The opposite limit corresponds to the situation that 
the two bonds with larger conductances in each triad 
dominate the charge transport. 

For a constant DOS, one obtains straightforwardly the 
number n{ei) and then the probability Pn{^i)- Substi- 
tuting them into Eq. (42) we finally obtain {R^j'jr')c — 



so-osm. and 



e7^i'')c ^ e° "8^«- (see SI for details). The longitudi- 



Dmax 

Sjk /C — C ^ , \C,j^j. 

'^ijk /c 

nal conductivity is obtained based on the 2-site func- 
tion Gij which should be no less than Gc in a perco- 
lation path. The result of axx equals Gc divided by 
the correlation length of the network and takes the form 
<^xx — o'o(T)e~^^'=, where cro(T) gives at most a power- 
law on T in [25]. Comparing this form with the AHC, we 



reach {a^"} 



'la/b 7a/!, 



Uxx with 7a = 1.76 and 



7b = 1.38. This leads to the scaling relation, the central 
result of this Letter, between and axx of the AHE 
in the insulating regime: 

1.38 <7< 1.76. (4) 
The maximum (minimum) of the AHC corresponds to 
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the smaller (larger) power index 7;, (70). This scaling 
range can be confirmed with a numerical calculation of 



the Eq. ( 42 1 . Furthermore, a direct numerical study for 
the configuration integral ^ gives the scaling exponent 
7 ft! 1.62, which is consistent with our prediction of the 
lower and upper limits. 

So far in the calculation we have assumed a constant 
DOS. This approximation is applicable for the ferromag- 
netic system with strong exchange interaction between 
local magnetic moments and charge carriers (e.g. oxides, 
magnetites) and half metals in general. In this case we 
do not need to sum over spin-up and spin-down states 
which contribute oppositely to the AHE, and the previ- 
ous results are valid. 

However, when the Fermi energy crosses both spin-up 
and -down impurity states, a symmetric DOS with p(e) — 
p{—e) leads to zero AHC. This is because under the trans- 
formation e;^CT — >■ —ei^-a- {I = i,j,k), Gijk changes sign, 
while Gij is invariant. Thus the averaging for AHC over 
all spin states and on-site energies cancels [T7]- We re- 
lax the previous simplifying restriction by expanding the 
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FIG. 4: (Color online) Scaling relation between the AHC and 
longitudinal conductivity. The theoretical results are com- 
pared with the experimental observations. 



DOS by Pie) = ^ 



n n! de^ 



e", where |e| < and we con- 



sider po = Pi'^p) > 0. Substituting this expansion into the Advancement of Science. 



An) 



with the 1st and 2nd 

dpo 



and 



Eq. (I2j) yields <^ = ELo"-y' 
nonzero terms respectively proportional to 
We can similarly evaluate the lower and upper limits of 
as before. The first two nonzero terms in the expan- 

^'^^ ■ ■ - ^'^^UT)cTl-^^''<y2i'' and 
'-ec{T)al-''^''a2Z^\ The 



SlOn are XCTxy \min/max 



rain I max 



Qm2M- 



appearance of M is due to the summation over the spin- 
up and -down states. We have also employed the result 
(|e|)c = 0.112$c. The specific formulas oiaQ{T) and £,c{T) 
do not affect the qualitative scaling between a-^^ and 
a XX- For the Mott and ES hopping regimes, we have re- 



spectively = kBT{To/TY'^ and = kBT{To/T) 
with To the constant depending on the DOS PlI^ET]. 
Note that ai^J and a^J have different physical meanings. 
The term ai^J dominates when the DOS varies monoton- 
ically versus e. Furthermore, when the DOS has a local 
minimum at the Fermi level, which may be obtained due 
to particle-particle interaction (coulomb interaction), we 



,1/2 



have dp/dep =^ 0. Then the term a^y varnishes and ax 
dominates the AHE. The above results also indicate that 
the AHC may change sign when dpo/dep or d'^po/dep 
changes sign, which is consistent with the observation by 
Allen et al^. 

Fig. [4] shows our theoretical prediction is consistent 
with the experimental observations of the scaling relation 
in this regime, hence completing the understanding of the 
phase diagram of the AHE. 
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SUPPLEMENTARY INFORMATION FOR "SCALING OF THE ANOMALOUS HALL EFFECT IN THE 

INSULATING REGIME" 



HOPPING MATRIX 



In the case the magnetization is saturated and thus M ~ Mcz, we rewrite the Hamiltonian Hp in the diagonal 
basis of the exchange term and obtain 

Hp — 'y ^ f-igC-igCia — ^ ^ ^ia,j P'^ia^j (3-, (5) 

where e^Q — ei + MTaa- Below are two different examples. First, for the dilute Gai_j;Mnj;As, the matrix Uajp 
describes the hopping of the holes localized on the Mn impurities. Under the spherical approximation tiajfj can be 
obtained based on by a unitary rotation [/(R^j ) from the direction to the hopping direction i ^ j [T]. We thus have 
Uajis = [U''(Rij)tdiagU{Rij)]a(3 with tdiag = diag[i3/2 , ^1/2 , t- 1/2 , ^-3/2] representing the situation that the hopping 
direction is along the z axis. Another case is for the localized s-orbital electrons. In this case, the hopping is given by 
tij = U^'Rij)[iij{l + iv^j ■ S)]U{Tlij). Here = diag[ti/2, ^.1/2] and v,.j ^ | /I" (VV^(r) x dr') with V(r) including 
the ion and external potentials, the spin-orbit coupling coefficient a = h/{4m^c'^) and m the effective mass of the 
electron. 



CONFIGURATIONAL INTEGRALS 



The averaging of a m-site physical quantity F{ei, 62, tm] ^ii ^2, ^m) along critical percolation path/cluster is 
given by 



1 



dei / de2... / dem x 



J dVi2 J d^f23-- J dV„,_i^„p(ei) Pfe(ei)/5(e2) x 

k — 7li 

00 00 

Pfe(e2)---p(em) ^ Pfc(em)F(ei,e2,...,em;ri,f2,...,rv„), (6) 



where Pn{^i,£,c) = {n-iy. /o"*''^'^ ^ ^'^^ H- Some examples are given below. The first one is the average value of 

n(ei,^c) in the percolation cluster. Note n{ei,^c) is a 1-site function. The averaging is straightforward and 

/ dein{e^)p{ei)J2T=iPkiii) _ I dein{ei) p{ei)n(ei) 



/*p(e»)E^i^/c(ei) J de^n{^i)p{e,) 



(7) 



The hopping conduction occurs when the average value n reaches the critical value ric- When the DOS p{ei) — po is 

Po (f: _ 

2^ PO s{L-\u\nec~\^^??de 



a constant, the number n(ei) is given by n{ei) — ^ (2o^Tp"(^c ^ ki|)^('?c ~ kiP)- Then we have 



3 {2akBTY /(^,-le,l)2(42_|g^|2)rf,^ 
from which we obtain the cut-off value ^j, by 



(2akBTy 



(2akBTfn, 
0.4067rpo 



1/4 



(8) 



(9) 
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Thus it gives 



Kc^[^Y\ To = 16^, (10) 



which is the Mott law. Accordingly, if we assume the density of states p(e) ^ e^, we obtain straightforwardly the 

Efros-Shklovskii (E-S) law =(^)^^^ M- Second, we give the formula for the longitudinal resistance based on the 
2-site function Zij — l/Gij. The longitudinal resistance for a percolation path is calculated by 

where N is the number of links along the percolation path. The above formula can be simplified by the fact that 
Sfcl2 Pk{ei) = n{ei,^c) - fi(ei,Cc) oc [n{e,,£,c)f. We then reach 

^ ^ N J dej J dej J d^njZjjje,, ej;nj)piei)[n{e„ ^c)? p{ej)[n{ej , ^c)? ^^^2) 
J de,J dcj J d^fijp{ei)[n{e^,^cWp{ej)[n{ej,(aW 

The longitudinal resistivity is given by Rxx/in-dLx), with the density of the percolation paths and the length of 
the material along x direction [2^ . Finally, if the physical quantity is a function of a triad with each site of the triad 
having at least three sites connected to it, the averaging of such physical quantity is given by 

F{ei £2 eg-ri fa fg) = / c^eic^eadeg / d^fia / d'^fagf (e; f)p(£i)[n(£i)]3p(e2)[n(£2)]^p(£3) [^^(es)]^ ^^3^ 

/ d£irf£2d£3 / d3fi2 / d3f23p(£i)[n(£i)]3p(£2)[n(£2)]3/9(£3)[n(£3)]3 

The anomalous Hall conductivity/resistivity will be calculated with this formula. 



FORMULA FOR MACROSCOPIC ANOMALOUS HALL CONDUCTIVITY 

Now we show rigorously the formula for macroscopic AHC in the hopping regime. The transverse voltage difference 
for the region from x — Ax to a: + Ax (Fig. 3 in the manuscript) reads 

Vy{x)^V,'' + V," + ... + V^. (14) 

For the general situation we allow some V/^'s to be zero (see Fig. 3 in the manuscript). In that case no triad 
forms for the incoming current li under the condition all direct conductances in a triad must be no less than Gc- 
To calculate V/^ , the voltage contributed by the i-th triad, we employ perturbation theory to the equation [4] 
lij = GijVij + Gijk(yik + ^jk)- First, in the zeroth order, we consider only the normal current, namely, the Hall 
current is zero and thus J2j hj — J2j ^ij^iP — 0; with which one can determine the voltage l^-^^-* at each site. Then, 
for the first-order perturbation, we have J2j ^ij — + J2j ^if^ ^ 0' which leads to jf'^^ — J2j '^if^ — 

J2j J2k GijkiVjk^ + Vik^) = ^ J2j GijVij. The current can also be written as 

= E E ^^.^(^^' + y!^) = \T. + - - = 2 E s.^l^- d^) 

j k jk jk 

For the hopping regime, the triads are dilutedly distributed and the Hall voltages induced by different triads are 




FIG. 5: (Color online) Resistor network transformation. 
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considered to be uncorrelated. Therefore, we obtain the Hall voltage of the i-th triad from the transformation indicated 
in fig. [5] that 



^ZlZ2'-^*2*3 ^ ^^1*3^^2*3 ~^ ^23*1 '-^'tl 't2 ^'tl't2^*2*3 ^ ^?1?3^*2*3 ~^ ^?3ll^2l22 



From the resistor network configuration one can see ^« ~ ^/q- For convenience, we denote li — 2IoXi{x) with 

J2i = 1- Generally (x) — J2i ^i^^ is a function of position x, and one needs to average it along the x direction. 
For a macroscopic system, one has N{x) — > oo. Furthermore, we consider at the position x, for each Ai there are 
ni{x) number triads that have such same current fraction A.^. Thus we have 

-, p qU) 

= 6/0^ / E E G G +GG +G G ^ 

J , , ._, ^7172""12J3 ""JU3""J2J3 I ""J3jl"-^jlj2 



To simplify this formula we extend the current distribution {A^} for the region between x — Ax and x + Aa; to 
the whole space along x direction, and then we can exchange the order of the integral and the first summation: 
J dx'^^^,y -^iSJlT^ ~^ ^iJT I dx'^^l}!^\ In the hmit N{x) — >■ oo and the length much larger than 

the typical length L of the triad, the calculation J dx X]j=f gives the average of all possible configurations of the 
triads through the percolating cluster. This leads to 

1^ J ^ili2""«2i3 ~r "-^iiis '-'2223 ~r U"i3ii u-jj 22 

with Ui — (l/Lx) J dxni{x) the average number of triads with in/outgoing current li. Note the identity ^^UiXi — 1 
is independent of position x, and therefore we have also '^^niXi = 1. The transverse electric field is given by 
Ey = /Ly. The longitudinal current density reads jo = Io/(LyL), where LyL represents the area of the cross 
section. With these results we obtain the Hall conductivity 

= oLa. 



xy xx\ f-i . f~< , y~i f~i 

^Ill2'^l2l3 ~ ^2123^22*3 ' "-^2311^1112 



'^^'^xx-^{ ,2^(2)^(2) , ,, , |2r(2)T(2) , ly- + ,2T^2)r^(2)>'=^ ^^^) 



where T^^' and T^jl are defined by 



with Ajj = e^Q, — , and 



T^^-* = lA . . A |g-2fc^(l<^3/3| + l<^'=Tl + l<^'=-<^'=Tl + l<^'=-<^3/3|) 

iik I 2.7 ^k\ 

^ijAjfele 



ijk ~ I'-^ij'- 

A- -A |e^2fc5T(l'^i°l + l^'=Tl + l<^i/3-«fc^l + l«ia-«j/3|) 



+ |AifcAfcj|e"5^^''-°'+''^'''+'''°""=^'+'"=^"'^'''^ (21) 



The configuration integral will be performed according to the Eq. (13 1. 



UPPER AND LOWER LIMITS 



For the lower limit, I < Icj/sl < \ek-y\- By keeping only the maximum term in the 

denominator and the minimum one in the numerator of the Eq. ( 19 1 we obtain 

tmax 
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To make the calculation realistic, we further consider the approximation by replacing the configuration integral of the 
exponential functions by configuration integral of the exponents. Then we get 



Similarly, the upper limit can be formulated with the restrictions Rij,Rjk > Rik and \eia\ > {(-j/sl > |eA:7l- By the 
same procedure we obtain 

Lower limit 

First we calculate the lower limit of AHC, which is given by 

k T 

{^xy }mm — '^^^ xx ^ (g) (-^ijlT) ci^Tjk^) (25) 
C tmax 

where (i?^r)c = e'^<^-^+^^-'=-«"=>= U„.,fl,,<i?,, , (e™"), = e0-5^'<l^'l+l^.l+l^.-''=l-l^»-^''l>=||,,l<|,^l<|,^|. We neglect the 
spin indices. The configuration integral {Rij + Rjk — Rik)c\Rij.Rjk<Rik given by 

(R +R -R ) = I '^^^'^''j'^^k f d^Rtj J d^R3kp{ez)Hei}fp{€j)[n{€j)fp{^^^^ + Rjk - Rik) , . 

^ /de,de,defc/d3R,,/d3R^.,p(e^)[n(e,)]3p(e^.)[n(e,)]3p(efe)[n(efe)]3 ' ^ ' 

with Rij, Rjk < Rik- Wc shall first perform the integral over position / d^Rij J (PUjk- Let R^j = Ri, R^fe = R2, and 
then i?3 Rik = Rf + R2 — 2R1R2 cos 9. Denote by the integral I ~ ^ j (P'R.ij J d^Iljk{Rij + Rjk ~ Rik) with 
J\fr = J (Ptiij J (PRjk- To write down the explicit formula of this integral, we apply the restrictions: Ri < Ri^max 
and i?i,i?2 < ^3, with Ri.max determined through 2ai?^°^ + 5/3(|ei| + \ej\ + \ei — €j\) = /3^c (from the condition 
G^™ = Gc or Z^'J"'^ = 1/Gc). With the basic triangle geometry (Fig. |6| we obtain 

^^2 




(a = 71-20) 

FIG. 6: (Color online) Triangle geometry for the configuration integral over the position space. 



yStt^ [ dR2Rl[ [ dd [ dRiRf sin0{Ri + R2 - J Rj + Rj - 2R1R2 cos0) 

'•'r Jo Jtt/2 Jq V - 
I-7V/2 i-Rt I 

-/ del dRiRlsine{Ri+ R2- J Rl + RI-2R1R2 COB 6)], (27) 

Jit/?, J2R9COS0 



where 



Ra = inm{Ri„,ax,R2 + V A2 - i?2 sin^ 6*}, 

, 7?n COS 

2 COS 6* 



R I 

i?b = min{i?i,„„„;^-^,i?2Cos0+A/A2-i?2sin2 0}, (28) 



with 



A^ = ^[^-^(N + N + h-..|)]^ (29) 
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Therefore the integral domain is not uniquely specified and depends on the the integral variables, which makes the Eq. 
(27) be still not analytically solvable. We need to simplify it by amplifying the integral domain. From the geometry 
of the triangle composed of {Ri, R2, R3), we can show the following inequalities: 



^R, cos 6+^K^-Rlsin^e ^^^^2 ^.^ ^ 



> 



(30) 



which is needed in the case R2 cos 9 + y — R2 sin 6 < Rimax , and 

jvRi^ax ^^^^2 sin0(i?i + i?2 - v/i?f + - 2R1R2 cos 9) 
JI^^'""''' dRiRlsine 

& cos e '^^i^i sin ^(^1 +R2- ^jRl + Rl-2RiR2 cos 9) 
fil" sin 61 

'f2R2 cos ^ ^ i 

when i?b < Rimax- Based on these results, we find that 



(31) 



< tt^Stt^ / rfi?2i?2 [ / d9 [ dRiRf sin6l(i?i + i?2 - A/i?! + R2 - 2R1R2 cos 61) 

■'Vr Jo Jit/2 Jo * 

+ / d9 dRiRlsm9{Ri+ R2- J Rl + Rl-2RiR2Cos9)\, (32) 

J-nIi J2i?9Cose 



with 

rR2max pR\max pR2max /''^/2 rRlmax 

Afr = Sn'^[ dR2Rl d9 dRiRlsm9 + dR2Rl d9 dRiRlsm9]. (33) 

Jo Jit/2 Jo Jo Jtt/S J2_R2Cose 



Employing the integral J^j^ d9 sin 6i^i?f + i?i - 2i?ii?2 cos 6I) [(i?i + i?2)^ - (i?f + i?i - i?ii?2)^/^] , we get 

finally 

^ - " - 0.424^2i?Lx/A/; (34) 

with Rmax — ^SLx{Rimax , R27nax} ■ It is casy to obtaiu the normalization factor as JVr ~ jgT^'^Rmax- After the integral 
over position given above we can now do it over the on-site energies. This gives 

,p ,P p\| 0.424 /de,dejdefcp(e,)[n(e.)]^p(ej)[n(ej)]3p(efc)[n(efc)]^i?L. 

[Hi + it2 ~ -K3;c|_Ri,_R,2<-R3 



23/18 /de.de,defe/p(e,)[n(eO]V(ei)[n(e,)]3p(e,)[n(efc)]3i?^„^ 
~ 0.l56/3£,c/a. (35) 

In above calculation we have considered the approximation that the density of states is a constant. 

Now we evaluate the average of energy. Similarly, the configurational average of the energy is given by 

(l^^l + kjl + ~ - 1^^ - efc|)c||ei|<|ej|<|efc| = 

1 J de^dejd€kp{€^)[n{ei)fp{€j)[n{ej)fp{ek)[n{ek)f{\e^\ + |ej-| + [e^- - efc| - la - ek\) 
2kBT J deidejdekp{ei)[n{e,)]^p(ej)[n{ej)]^p{ek)[n{ekW 

(36) 

To simplify the above integral, we check \ej — ek\ — \ei — ek\ with the restriction: |ei| < \ej\ < \ek\- For the case i) 
Sgn(ej) = Sgn(ej) = Sgn(efc) = ±1, we have |ej -e^l - je,; -e^l = -\ei-ej\; For ii) Sgn(ej) = Sgn(ej) = -Sgn(efc) = ±1, 
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we have |ej-efc|-|ei-efc| = -\ei-ej\; For iii) Sgn(ei) = Sgn(efe) = -Sgn(ej) = ±1, we have |ej-efc|-|ei-efe| = -\e. -ej\; 
For iv) Sgn(ej) — Sgn(efe) = — Sgn(ei) = ±1, we have \ej — ek\ — \ti — = — Cjl- For this we obtain that 

+ h\ + h ~ efcl - V^ - efcDc - + kjl ^i\)c- (37) 

Then by a straightforward calculation one can verify that 

f + + - efel - V^ - efc|)c||.,|<|.,|<|.,| = 0.086/3Cc. (38) 



From eqs. (351 and (38) we have 



a(i?,y - i?.fc)U,,,fl,,<i?,, + ^(|e.| + Vj\ + - efcl - |e, - efe|)c||e.|<|.,|<|..| = 0.242/3^^. (39) 

The lower limit of the AH conductivity is then obtained by 

X^xy imin = 3La^^ ^ 6 (40) 



The longitudinal conductivity is given by Gc divided by the correlation length of the network and thus takes the 

To' 

T , 



form = o'o(2^)e (for the Mott hopping regime, one has /3^c =( t^)^^^)- We reach further 



jc^xy Imm 3Lcr(, ^ (q) ^xx «^xx> 7-1-76. (41) 

6 67; 



^raax 



UPPER LIMIT 



Now we show the result of the upper limit, which can be done in a similar procedure. The upper limit is given by 

AH 2 ^bT 

Wxy }max — Sicr^.^ {Rijk^} ci^ijk^} c, (42) 

where (i?™r)c = e"<'^*^+'^-^'="'^'""^^ |fl.„fl,.>fl.. , (e^Dc = e"-S^<l'-l+l'^l+l^-^'^-l-l^'-^'=l)=||e.|>|e,|>|.,|. To calcu- 
late {Rij + Rjk - Rik)c\Rij,R,k>Rik we again consider first the integral I ^ d^R^ J d^R2{Ri + R2 - 
y^Rl + R2 — 2R1R2COS9) with Afr = J (PRi J d^R2. Note the integral restrictions for the upper limit are: 
Ri < Ri,max and i?i , i?2 > ^^3, and with the triangle geometry (Fig. ([6])) we obtain 

/ = ^Stt^ / dR2Rl j de I dRiRl sin6l(i?i + R2 - Rj + Rl - 2R1R2 cosO), 

(43) 



2 003 e 



where Ra = mm{Rimax,2R2Cos6,R2COs9 + yA^ — i?|sin^0}. Again we simplify the integral by amplifying the 
integral domain. For this we consider the following inequality: 

gcTsVyx^^i^p^ dRiRl sin ojRi + R2~ VRj + ^RiR2 cos e) ^ 

jR. cos e+VA^-flisin^fl ^^^^2 ^(^^ ^ _ + i?2 _ 2i?ii?2 COS 0) 



(44) 



which is needed in the case R2 cos 9 + JA^ - Rl sin^ e < Rimax- With this we find that 



I-j^^ir'^j dR2Rl j de J ^ dR^Rl sm d {Ri + R2 - ^ Rl + R:1- 2R1R2 con 9). (45) 
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By the same procedure used in the lower Uniit we obtain / ~ 0.37297r^i?J„„,j./A/'r, and Mr = 
^tt"^ dR2R\ J^^^ dO dRiR\s\ii9 = 0.36l7r^i?^^^. Further doing the integral over the on-site energies 

yields {R^j + Rjk - Rik) c\rI^X>^<r,, = QASZKc/a.. 

The configurational average of energy + + — efc| — — efe|)c||ei|>|e |>|efc| can be simplified by checking 

jci — efel — \ej — efcl with the restriction: \ei\ > \ej\ > \ek\. Through a similar analysis as applied in the lower limit one 
can verify ( | Cj | + | | + | — | — | Cj — | ) c ~ ( | ei | + | | + 1 1 e j — ej | ) c ■ Substituting this result into the original integral 
we obtain finally i^(|ei| + |e,| + |e,- - efc| - \ei - efe|)e||e,|>|e,|>|e.| = 0.1375/3^^- For this we obtain (i?^r)c ^ e0-483^«s 
{e^r)c - e0-i3^^«<=, and the upper limit of the AHC by 

S^^Hx or ^0.621 1.379 ^ „, ^ 1 on /A(y\ 

Wj; Imax = 3Lcro ,q. tT^^ OC CT^'^, 7 ~ 1.38. (46) 

Based on the results obtained above we thus conclude {cr^^} oc a].^ with 1.38 < 7 < 1.76. 
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